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Abstract. In this paper, we construct stationary classical solutions of the incompress- 
ible Euler equation approximating singular stationary solutions of this equation. This 
procedure is carried out by constructing solutions to the following elliptic problem 

f -e*Au = J2T =1 Xn, (u - q - % In §)* , x e Q, 
\u = 0, x€dQ, 

where p > 1, ft C M 2 is a bounded domain, fli CC fl, i — 1 • • • , m are suitable small 
domains such that £li f] £lj is empty if i ^ j and q is a harmonic function. 

We showed that if fl is simply-connected smooth domain, then for any given stable 
critical point of Kirchhoff-Routh function W(zi , • • • ,z m ) with the positive strength Ki > 0, 
there is a stationary classical solution approximating stationary m points vortex solution 
of incompressible Euler equations with vorticity Y^JLi K i- 

Existence and asymptotic behavior of single point non- vanishing vortex solutions were 
studied by D. Smets and J. Van Schaftingen in |22) . 
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1. Introduction and main results 
The incompressible Euler equations 

\ t + (v • V)v = - VP, 

V-v = 0, 1 ' ' 

describe the evolution of the velocity v and the pressure P in an incompressible flow. In 
1R 2 , the vorticity of the flow is defined by u> = V x v := diV% — c^fi, which satisfies the 
equation 

uj t + v • Vcj = 0. 

Suppose that u is known, then the velocity v can be recovered by Biot-Savart law as 
following: 

1 -x 1 - 



V = u * 



2n \x\ 
l 
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where x 1 - = (x 2 , —x\) if x = (xi, x 2 ). One special singular solutions of Euler equations is 
given by u = YlT=i K ^Xi(t), which is related 

EKi {x — Xj(t)) -1 " 
, 2^ \x-xAt)\ 2 ' 
i=i 1 v ' 1 

and the positions of the vortices X{ : R — > M 2 satisfy the following Kirchhoff law: 
where W is the so called Kirchhoff-Routh function defined by 

m 

2 27T SCj — X 



:i l 



In simply-connected bounded domain ft C M 2 , similar singular solutions also exist. Sup- 
pose that the normal component of v vanishes on <9ft, then the Kirchhoff-Routh function 
is 

^ m ^ m 

■ • • , SC m ) = - 2j KiKjGfa, Xj) + - 22 f^H(Xi, Xi), (1.2) 

where G is the Green function of —A on ft with Dirichlet boundary condition and H is 
its regular part (the Robin function). 

Let v n be the outward component of the velocity v on the boundary <9ft, then we see 
that J dn v n = due to the fact that V ■ v = 0. Suppose that v is the unique harmonic 
field whose normal component on the boundary <9ft is v n , then vo satisfies 

V • v = 0, in ft, 

V x v = 0, in ft, (1.3) 
n ■ v = v n , on 9ft. 

If ft is simply-connected, then v can be written v = (Vil ! o)~ L , where the stream function 
ipo is determined up to a constant by 

-Aipo = 0, in ft, 
v n , on o\l, 



Or 

where denotes the tangential derivative on <9ft. The Kirchhoff-Routh function associ- 
ated to the vortex dynamics becomes(see Lin [17] ) 



W(xi,--- ,X m ) = ^ ^ KjKjGjXj, Xj) + n ^ K?H(Xj, Xj) + ^ Kj^Xj) ■ (1.5) 
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It is known that critical points of the Kirchhoff-Routh function W give rise to stationary 
vortex points solutions of the Euler equations. As for the existence of critical points of W 
given by (II. 2ft . we refer to j5]. 

Roughly speaking, there are two methods to construct stationary solutions of the Euler 
equation, which are the vorticity method and the stream-function method. The vorticity 
method was first established by Arnold and Khesin [3] and further developed by Burton 
[7] and Turkington [23] . 

The stream-function method consists in observing that if ip satisfies — Aip = f(ip) for 
some function / E C 1 (M), then v = (VV ; ) ± and P = F(ip) — ^\Vip\ 2 is a stationary 
solution to the Euler equations, where (V?/')~ L := (Jj- ; — Jr")'-^^) = f f(s)ds. Moreover, 
the velocity v is irrotational on the set where f(ip) = 0. 

Set q = —ipo and u = ip ~ i>o, then u satisfies the following boundary value problem 



In addition, if we suppose that inf^ q > and f(t) = 0, t < 0, the vorticity set {x : f{jp) > 
0} is bounded away from the boundary. 

The motivation to study (II. 6p is to justify the weak formulation for point vortex solu- 
tions of the incompressible Euler equations by approximating these solutions with classical 
solutions. 

Marchioro and Pulvirenti [IH] have approximated these solutions on finite time intervals 
by considering regularized initial data for the vorticity. On the other hand, the station- 
ary point vortex solutions can also be approximated by stationary classical solutions. See 
e.g. [3 [211 [22], [231 121] and the references therein. It is worth pointing out that the above 
approximations can just give explanation for the formulation to single point vortex solu- 
tions. In this paper, we will show that multi-point vortex solutions can be approximated 
by stationary classical solutions. There are many results for problem (11.61) on the existence 
and asymptotic behavior of solutions under various assumptions. In (H El [121 [2D [23], the 
constrained variation methods were used to find solutions for the equation 



under the constraint f n F{u — q) = fi, where A > is a Lagrange multiplier a priori 
unknown. On the other hand, in [H [21 [201 [2H [2S], the solutions were obtained by using 
Mountain Pass Lemma for various nonlinearities. For the asymptotic behavior, Berger and 
Fraenkel [6] began studying the asymptotic behavior for variable fi and q, but the lack of 
information about A is still an obstacle. 

To avoid this obstacle, Yang [21] studied the minimization of the functional | J n \ Vu\ 2 — 
\ f n F(u — q) under the natural constraint J n |Vw| 2 — ^ J n uf(u — q) = and obtained 
the asymptotic behavior of the solutions u £ as e — > for Q = R^, q(x) = Wx\ + d, where 
W } d > 0. That is, set A £ = {x E M+ : f(u £ —q) > 0}, k £ = J n f{u £ -q) and x £ E A £ , then 




Au = f(u - q), 



x E f2, 
x E dQ. 



(1.6) 




x E d£l, 



(1.7) 
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diamA £ -> 0, dist{x e ,dR%) and j£ - G(x e , •) ->■ in W/ ' c r (M^) for r G [1,2). Later 
on, similar results were obtained in [T5] for bounded domains with additional information 
that q(x £ ) — > min$-2 q. However, it has been pointed out in |22J that the solutions obtained 
above corresponded to desingularization of point- vortex solutions with vanishing vorticity. 
To get non-vanishing vortex solutions, D. Smets and J. Van Schaftingen [22J investigated 
the following problem 

-e 2 Au= (u-q-%\n\) p + , in fi, 
w = 0, on <9fi, 

and gave exact asymptotic behavior and expansion of the least energy solution by esti- 
mating the upper bounds on the energy. The solutions for (11.91) in [22] were obtained by 
finding a minimizer of the corresponding functional in a suitable function space, which can 
only give approximation to a single point non-vanishing vortex. This method is hard to 
obtain multiple non- vanishing solutions. 

In this paper, we approximate stationary vortex solutions of Euler equations fll.lj) with 
multiple non-vanishing vorticity by stationary classical solutions. 

Our main result concerning (11. ip is the following: 

Theorem 1.1. Suppose that Q C R 2 is a bounded simply- connected smooth domain. Let 
v n : dfl — > R be such that v n G L s (dQ) for some s > 1 satisfying j dn v n = 0. Let 
Hi > 0, % = 1, • • • , m. Then, for any given stable critical point (x*, • • • , x£j of Kirchhoff- 
Routh function W(xi, • • • ,x m ) defined by (ll.5p . there exists e > 0, such that for each 
e G (0,£o); problem (II. ip has a stationary solution v e with outward boundary flux given by 
v n , such that its vorticities u e satisfying 

supp(u £ ) C Uf =l B(x ii£ ,Ce) forx it£ eQ, i = 1, • • • ,m 

and as e — > 




Remark 1.2. The simplest case, corresponding to a single point vortex ( m — 1 ) was studied 
by Smets and Van Schaftingen [22] by minimizing the corresponding energy functional. 
In their paper W(xi )£ ) — > sup xe nW(x). Even in the case m — 1, our result extends 
theirs to general critical points (with additional assumption that the critical point is non- 
degenerate). The method used in [22] can not be applied to deal with general critical point 
cases. The method used here is constructive and is completely different from theirs. 

Remark 1.3. In this case that m — 1 suppose that x% is a strict local maximum(or mini- 
mum) point of Kirchhoff-Routh function W(x) defined by (jl.5p . statement of Theorem ll.il 
still holds which can be proved similarly(see Remark 1 1.51) . Thus we can obtain correspond- 
ing existence result in [22J. 



REGULARIZATION OF POINT VORTICES 



5 



Theorem 11.11 is proved via considering the following problem 



f-5 2 A M = E7 =1 x^(«-g-iini) 

\u = 0, 



(1.9) 



where p > 1, q G C 2 (Q), Q C M 2 is a bounded domain, fij C is a subdomain such that 
x* £ Qj, j = 1, ■ ■ ■ ,m and fij fl Qj — if i ^ j. 

Theorem 1.4. Suppose q G C 2 (Q). For K{ > 0, i = 1, • • • ,m. Then, for any given 
C 1 -stable critical point (x*, ■ ■ • , x^J of Kirchhoff-Routh function W(xi, ■ ■ ■ , x m ) defined by 
fll.5p . t/iere exists £ > 0, such that for each e G (0,e ); dE5]) a solution u £ , such that 
the set fl e>i = {x : u £ (x) — ^- In i — q(x) > 0} CC Of, i = 1, 2, • • • , m anc? eac/i fi £j j shrinks 
to x* eQ, as e ->• 0. 

Remark 1.5. For the case m = 1 suppose that xi is a strict local maximum(or minimum) 
point of Kirchhoff-Routh function W(x) defined by (II. 5ft . statement of Theorem 11.41 still 
holds which can be proved by making corresponding modification of the proof of Theorem 
11.41 in obtaining critical point of K(z) defined by (I4.ip (see Propositions 2.3,2.5 and 2.6 [9] 
for detailed arguments). 

For domains which may not be simply-connected, we show in the following result that 
the topology of the domain plays an important role in the existence of solutions. 

Theorem 1.6. Suppose that the homology ofQ is nontrivial. Then, for any positive integer 
m, there exists Eo > 0, such that for each e G (0, Eq), (11.91) has a solution u £ , such that the 
set Q £t i = {x : u e (x) — |j In | — q(x) > 0} CC fij, % — 1,2, • • • ,m and each shrinks to 
a point x* G Q, as e — > 0. Moreover x* ^ x*, ifi^j. 

Remark 1.7. Since m is arbitrary, from Theorem II. 6| we can see that the number of 
solutions for (11.101) is unbounded as 5 — > 0. 

Not as in [22] where (ll.9p is investigated directly, we prove Theorem II. 1[ Theorem 11.41 
and Theorem 11.61 by considering an equivalent problem of (II. 9p instead. Let 

and 5 = e(jY !L j) S ^~ , then (II. 9ft becomes 



We will use a reduction argument to prove Theorem 11.41 and Theorem 11.61 To this end, 
we need to construct an approximate solution for (11.101) . For the problem studied in this 
paper, the corresponding "limit" problem in M 2 has no bounded nontrivial solution. So, 
we will follow the method in [TUJ [TTJ to construct an approximate solution. Since there 
are two parameters 5, e in problem (ll.lOp . which causes some difficulty, we must take 
this influence into careful consideration and give delicate estimates in order to perform the 
reduction argument. For example we need to consider (sij, • • • , s m ,<5) and (ai,,5, • • • , a mt s) 
together in Lemma 12.11 




(1.10) 



6 DAOMIN CAO, ZHONGYUAN LIU, AND JUNCHENG WEI 

As a final remark, we point out that problem (11.101) can be considered as a free boundary 
problem. Similar problems have been studied extensively. The reader can refer to [3, QUI 
[TTl [T3| H~6] for more results on this kind of problem. 

This paper is organized as follows. In section 2, we construct the approximate solution 
for fll.lOp . We will carry out a reduction argument in section 3 and the main results will 
be proved in section 4. We put some basic estimates in the appendix. 

2. Approximate solutions 
In the section, we will construct approximate solutions for fll.lOp . 

Let R > be a large constant, such that for any x G Q, Q CC Br(x). Consider the 
following problem: 

' -6 2 Aw = (w-af + , mB R (0), 

w = 0, ondB R (0), [ ' ' 

where a > is a constant. Then, (12. ip has a unique solution Ws, a , which can be written as 

M ^ \aln§/lnf , s 5 < |ar| < R, K } 

where <f>(x) = <fr(\x\) is the unique solution of 

and ss G (0, i?) satisfies 
which implies 



5| ln<5|0-i)/: 

Moreover, by Pohozaev identity, we can get that 



J Bi(0) Z JB^ 



p = 2tt|0'(1)|. 

i(0) * JBi(0) 

For any z£fl, define Vl^^x) = W^aC^ ~ z)- Because Ws )Z>a does not vanish on dfl, we 
need to make a projection. Let PWs tZ ,a be the solution of 

-5 2 Aw = (W 6 , z ,a - a) p + , in Q, 
w = 0, on <9f2. 



Then 



where g(x, z) satisfies 



PW StX ,a = W 5 , z , a - ^r9{x, Z), (2.3) 
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-Ag = 0, in O, 

R 
\x-z\ ' 



g = In T-r^-r, on d£l. 



It is easy to see that 

g(x, z) — In R + 2nh( y x, z), 

where h(x,z) = —H(x,z). 

We will construct solutions for fll.lOp of the form 

m 
0=1 

where Zj G Q for j = 1, • • • , m, oj$ is a perturbation term. To obtain a good estimate for 
W(5, we need to choose a$ 7 j properly. 
By (12. 3p . we have 

m m / m „ \ P 

- 5 2 A )T P^,,, aij - £ Xn, £ /,U — - " " 77 ' 7 



j=l j=l \i=l 



In el 



£ (WW,,, - a Slj ) p + - J2 xn 3 A Yl ( w &wi ~ y^kaiv, *i) I - «i 



j=i J=l \i=l 



ln ii — i ' lne 



(2.4) 

Denote Z = (z\, • • • , z m ) G R 2m . In this paper, we always assume that Zj G £1 satisfies 

d(zj,dti) > g > 0, - > £ L , i, j = 1, • • • ,m, i ^ j, (2.5) 
where g > is a fixed small constant and L > is a fixed large constant. 

Lemma 2.1. For S > small, there exist (ss,i(Z), ■ ■ ■ , s<5 >m (Z)) and (0,5,1 (Z), • • • , a,5, m (Z)) 
satisfying the following system 

s ^- % -^ m = _^_ y i = ( , 6) 

and 

a < = Ki + Thi^ + ln fl a * " ln fl a i' i = 1, • • • , m, (2.7) 



where for % = 1, • • • , m ; 



G(x, z,-) = In I — - — - - g(x, Zi). 

\X Zj\ 
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Proof. We will show that system (I2.6p - fl2.7p has a solution (si, • • ■ , s m , Oi, • • • , a m ) in V = 
[r^, 8\ ln5\ } m x n™Jf , It is easy to see, for fixed (s x , ■■■ ,s m ) with < 5 < 5* small, 

that (12. 7p has a solution (cti, • • • , a m ) depending on (si, • • • , s m ), such that y < < 
For such ( ) and for % — 1, • • • , m define 



2 

/ / / -i \ 

2 



0iOi, • • • , s m ) = -^-p h 



then it is easy to verify that 

0i(si, • • • ,s m ) > 0, Si = S\ \aS\,Sj e [-J-tt, 5\ ln<5| ] for j = 1, • ■ • ,m,j 7^ z 



ln<5| ' 

0i(si,--- ,s m ) < 0, Sj = Sj <E [p^j, <y|ln<J|] for j = 1, • • • ,m,j ^ i. 

By the Poincare-Miranda Theorem in [HI [TS], we can get (s^i, • • • , s«s, m ) such that 
0i( s <5,i; " • • > s 5,m) = 0. Therefore we have completed our proof of Lemma I2TT1 

□ 

For simplicity, for given Z = {z\, • • • ,z m ), in this paper, we will use as t i,sg t i instead of 
as,i{Z),s s>i (Z). 

Remark 2.2. More precisely, we have the following relation 

1 1 



In i n « 

ss,i £ 



In 


lne 


Ine 


2 



TibET ,»=l,--,m, (2.8) 



dz jjh V I m e l / ^i,h V I ln e \ J 

Indeed, ( 12. 8 j) can be deduced from (I2.6p (see [IT] , for example). (I2.9P can been deduced 
from (I2.8P and (12. 7p . Differentiating both sides of ( 12. 6 p and (12. 7p with respect to z^ h we 
can get a linear system of and jr^, which will deduces ( 12.1 OR . 

From now on we will always choose (05,1, • ■ ■ , a$, m ) and (s^i, ■ ■ • , s<5, m ) such that ( 12. 6 p 
and (12. 7p are satisfied. For (a^i, • ■ • , a^ m ) and (s^i, • • • , s<5, m ) chosen in such a way define 

Ps,z,j = PW s , Zj , as> . , P 5 ,z = Ps > z >r ( 2 - U ) 

3=1 

Then, we find that for x e BL Sgi (zi), where L > is any fixed constant, 
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Ps,Z,i( x ) ~ K i ~ 2 ^ X ) = W S , Zi ,a s>i (x) - ^h~9( x , z i) ~ «i - 



2nq(zi) 2tt 



lne| 



ln^ 



a S,; 



ln-^ 



(Dg(zi,Zi),x- + 0(\x - Zi\ 2 ) 



|lne| |lne| 



((Dq{zi),x - Zi) + 0{\x - Zl \ 2 )) 



2nq(zi) 2tt 



|lne| |lne | 



(Dq(zi),x - Zi) 



as,', 



as. 



In 



R 



g(z i: - - — ' w (Dg(z i} z { ),x- z { ) + O 



In -5 



In el 



and for j ^ i and a; G B Lsgi (zi), by @2J2|) 



In 

So, by using ( 12. 7ft . we obtain 

27rg(x) 



4,i \ 

|lne|A 



-P<5„z(a0 - Ki - 



|lne| 



=W5 lZi , ai .(a;) - a s ,i - t-^-t (Dq(zi),x - Zi) - 7^%- (Dg(z h z { ), x - ^) 



lne| 



In 



S5,; 



+ J2^(DG(zi,z j ),x-z i ) + 

3+i s S,j 



In e| 



x e B Ls (z^ 



(2.12) 



We end this section by giving the following formula which can be obtained by direct 
computation and will be used in the next two sections. 



dW 5>Zi , as Ax) 



dz ith 
< 1 



as.i \ (p +1 )/ 2 ,i/\x — Zi\\ Zih — Xh 
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3. THE REDUCTION 

Let 



W (X) 



<t>{\x\), \x\<l, 
In \x\, \x\ > 1. 



Then w G C 1 (IR 2 ). Since <f)'(l) < and In \x\ is harmonic for \x\ > 1, we see that w satisfies 

-Aw = w p + , inM 2 . (3.1) 
Moreover, since w + is Lip- continuous, by the Schauder estimate, w G C 2 ' a for any a G 
(0,1). 

Consider the following problem: 

- Av - pw p +\ = 0, v G L°°(M 2 ), (3.2) 
It is easy to see that i — 1,2, is a solution of (13. 2p . Moreover, from Dancer and Yan 
we know that w is also non-degenerate, in the sense that the kernel of the operator 
Lv := -Av-pw^v, v G D 1 ' 2 (R 2 ) is spanned by {§^ } ||}. 
Recall that Z = (z±, ■ ■ ■ , z m ), and Zj G £1 satisfies 

d(zj, dil) > q > 0, \zi — Zj\ > g L , i ^ j, (3.3) 

where g > is a fixed small constant, and L > is a large constant. 
Let Ps,z,j be the function defined in (12.111) . Set 



F s>z = { u : u G L P (Q), I = 0, j = 1, • ■ ■ , m, h = l,2 

In c,z j,h 



and 



E 5 ,z = : u G W 2 ' P (Q) n H^Q), jf A 

For any w G £ p (fi), define Q^w as follows: 

m 2 



u = 0, j = 1, • • • , m, h = 1, 2 



*"=-£E^(-^(^)) 



3=1 h=l 

where the constants bj^, j — 1, • • ■ , m, h — 1, 2, satisfy 



Since J n ^^-Qsu = 0, the operator Q,5 can be regarded as a projection from L p (Vt) to 
F^z- In order to show that we can solve (13.41) to obtain bj^, we just need the following 
estimate ( by fl2TT0|) and flgZEgp ): 
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2 f A ( dPs,z,j \ dP 6 , Zii 



n 



^ijhfc- r— r + O 



hi£\p +1 \\\ne\P +1 



e 



where c > is a constant, = 1, if i = j and h = h; otherwise, = 0. 
Set 



(3.5) 



m ( 27rg(x)V 
L 5 u = -5 2 Au - ^PXn, \P&,z - Kj |hi7j/ 



p-i 

u. 

3=1 v 7 + 



We have the following lemma. 



Lemma 3.1. There are constants p > and 5 > 0, such that for any S e (0, S ], Z 
satisfying (I3.3P . u G rni/i Q5L5U = m 0\ U r f L =1 B Ls5 ^ (zj) for some L > /arge, t/ien 



||Q<5 -^(5^11^(0) > ^° (p _ 1)2 |MU°°ffi)- 

I In 5 1 p 

Proof. Set s raj j = s$ n)3 -. We will use || • || • ||oo to denote || • \ and || • || respectively. 

We argue by contradiction. Suppose that there are 5 n — > 0, Z n satisfying ( 13. 3 p and 
«n e £a n ,z n with Q Sn L Sn u n = in f2 \ Uf =1 B LSn i:j {z j>n ), ||«n||oa = 1, such that 

2 

1 5n 

|lnd n | p 

Firstly, we estimate ^ n in the following formula: 



m 2 / dP' 

Qs n L 5n u n = L Sn u n - b i> h > n ( ~ 5 ™ A 

i=i fe=i ^ 



3- 

For each fixed z, multiplying (13 .6p by 8/ ^ n,z '"' , noting that 



6 n ,Z n ,j 

dzjh 



(3.6) 



f(Qs n L Sn u n )^f^ = 0, 
Jn ° z i,h 



we obtain 
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9Ps n ,z„,i\ _ [ [ T „, \ d Ps n ,z n ,i 



m 2 „ r 

3=1 h=l JQ V 



' 2 " ' dP, ,,\ dP. 



dz jh J dz ih 



Using (I2.12p and Lemma fA.lt we obtain 

dPs n ,z n ,i 



n 



'dP Sn>Zn A ^ f D , 2nq(x)\ p ~ 1 dP Sn ,z n 



f f A 2 a fdPs^A sr f P 27r g (x) V 



f ( W \P~1 / g ^»,Zi,n,af w .i _ 



dZi 



1 1 



n 



S n ,j \ \ O^S n ,Z n ,i 



j = l 3 

o 



dz itJl 



c 2 



Using (13. 5p . we find that 

h,h,n = (41 m£ n|) 

Therefore, 



i=i h=i v 



X 7, n 



dz 



j 

m 2 



P W {W Sn>ZjM - a Snjj ) + ^ — — j 



3=1 h=l 

£n \bj,h, 



m 2 ~1, 



°IEE 



it | 



1 tttS^ I "'(«)■ 



Thus, we obtain 
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Lx u„ = Qx Lx u„ + O I i — '—, — - |=0| —, — 5 

" " ! ln.:„ ; J \ n \lnS |^ 

For any fixed i, define 

Ui, n (y) = U n (s n ^y + Zi j7l ). 

Let 

f a sr <* (p ( + \ y- 1 

L n u = -Au - ^P-fitXQi \ P& n ,zS s n,iy + Zi,n) - «i - | lng i J 



u. 



Then 



Noting that 



we find that 



As a result, 



I si ~ 

S n,i * — 2~ \\I-'nU'i,n\\p = 1 1 L$ n U n \ \ p . 

S n,i 



8, ^ 2 



=0 



Sri. A J V I in (5^7 1^ 

2 

T . On 

Lxu n = o 



1 A 



L n Ui tTl =o(i), mL p (ft n ), 

where Q n = {y : s n<i y + z i;n G Q}. 

Since ||wi,n||oo — 1> by the regularity theory of elliptic equations, we may assume 

u^n-fui, inC/ oc (R 2 ). 

It is easy to see that 

<ST<i ( P ( + x 2vrg V- 1 

-p-xui yn n .zA s n,iy + Zi,n) - hi - | lng i J 



c 2 / / „ \\P-1 



J2 . W<5 n ,^,a^ - a 5nii + O I — ^ ) ) + o(l) 



p-i 
— >■ u>l . 



Then, by Lemma fA.lt we find that Ui satisfies 

— Aui — pw±~ U{ = 0. 
Now from the Proposition 3.1 in [IT] , we have 
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dw dw 

u i = c 1 - hc 2i -. (3.7) 

ax\ 0x2 



Since 



we find that 



y A( -a*r- K=0, 



rp-1 n 
?+ T—Ui = 0, 

OZh 

which, together with (13. 7p . gives Ui = 0. Thus, 

0, inC^S^O)), 

for any L > 0, which implies that u n = o(l) on dBL Sni (z itn ). 
By assumption, 

Qs n L 5n u n = 0, in O \ U t fe =1 S Xsn] .(zi,„). 
On the other hand, by Lemma IA.1| for j = 1, • • • , m, we have 



Ps n , Zn -^-^l) =0, xen,\B^(z j ,„). 



Thus, we find that 



-Au n = 0, mn\uT =1 B LSn Xz t ,n). 
However, u n = on <9f2 and u n = o(l) on dB LSnt (z itn ), i = 1, • • • , m. So we have 

w n = o(l). 

This is a contradiction. 

□ 

Proposition 3.2. Q^L,; is one to one and onto from E$ t z to F$ } z- 

Proof. Suppose that Q5L5U = 0. Then, by Lemma I3"TT| u = 0. Thus, QsL$ is one to one. 
Next, we prove that Q$L$ is an onto map from E$,z to Fg t z- 
Denote 

E= [u:ueH%(Sl), J D d Q*' Zd Du = 0, j = 1, • • • , m, /i = 1, 2}. 

Note that E s , z = E n H/ 2 ' p (fi). 

For any h G i^z, by the Riesz representation theorem, there is a unique u G Hq($Y), 
such that 



REGULARIZATION OF POINT VORTICES 



15 



5 2 [ DuDlp = f fop, V <p e H£(Q). 
Jn Jn 



(3i 



On the other hand, from h G F$,z, we find that u G E. Moreover, by the L p -estimate, we 
deduce that u G W 2 ' P (Q). As a result, u G E S<Z - Thus, we see that Qs{—S 2 A) = —5 2 A is 
an one to one and onto map from Eg,z to Fg t z- On the other hand, QgLgu = h is equivalent 
to 



u = P S~ 2 (-Q s A)- 1 



Qs Yl Xn > ( p ^ z 



d=i 



2irq(x) 
I In el 



p-i 



■ft 



+r 2 (-g 5 A)- 1 ^, « g £ 5)Z 



It is easy to check that 5 2 {—QgA) 1 
pact operator in Egz- By the Fredho 



Qs ( Y^j=\ Xn 3 (Ps,z 



3 |lne 

m alternative, ( 13. 9p is solvable if and only if 



(3.9) 
is a com- 



u = p<5- 2 (-g f A)- 1 



Qs \ Y XQ i ( Ps > z 
3=1 



2iTq(x) 
I In el 



p-i 



u 



has trivial solution, which is true since QgLg is a one to one map. Thus the result follows. 



Now consider the equation 



where 



Q5L5U = Qsh + QsRs(u 



h = 2^ Xn, [Ps,z ~ Kj - -— — I 
j—i \ 11/ 



( w s**u - a sj) + 

3=1 



□ 



(3.10) 



(3.11) 



and 



*«-) -<*+«- -&», - * - 

J 1 ~\ 5 1 + (3.12) 

/ P 2^) y- 1 

Using Proposition 13.21 we can rewrite (13 . 1 0[) as 

u> = G s co =: {Q 5 L 5 )- 1 Q 6 (k + Rs(u})). (3.13) 

The next Proposition enables us to reduce the problem of finding a solution for (ll.lOp 
to a finite dimensional problem. 
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Proposition 3.3. There is an 5 > 0, such that for any S G (0, S ] and Z satisfying (I3.3p , 
( 13.1 Op has a unique solution u$ G E^ z, with 

\\u3 S \\ 00 = o{s\\RS\^ S j. 
Proof. It follows from Lemma IA.ll that if L is large enough, 5 is small then 

( P6 ' z ~ Kj ~ Tm^) = °' x £V j \B Lss .(z j )J = l, 



Let 



M = E 5iZ n{||a;|| 0O <5|ln5| £ 2 i }. 



Then M is complete under L°° norm and Gs is a map from E$,z to ii^z- We will show 
that Gs is a contraction map from M to M. 
Step 1. G,5 is a map from M to M. 

For any u G M, similar to Lemma fA.lt it is easy to prove that for large L > 0, 5 small 



(p sz +u; - Kj - = 0, in \ B LsS j ( Zj ). 

Note also that for any u G L°°(Q), 

Q 5 u = u mn\VJj =l B LsS] {z 3 ). 
Therefore, using Lemma [A. 11 (13.111) and (13.121) . we find that for any u G M, 

Qsh + QsRsitu) = h + R s (u) = 0, in ft \ Uf =1 B Lssj { Zj ). 
So, we can apply Lemma [3.11 to obtain 

(p-i) 2 

\\{QsLs)- 1 {Q s ls + QsRsH)\\ 0O < ^^j " \\Qsl 5 + QsRs(u)\[ 
Thus, for any oj G M, we have 

Hoc ^(QsLs^Qsik + Rsi^Woo 



(3.14) 



(3.15) 



It follows from ( I3.4p ~ (j3.5p that the constant bj^, corresponding to u G L°°(ft), satisfies 



< C|ln5r +1 ^ [\ 
. T Jn l 



i, h 

Since 



\u\. 
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l s + R s (u) = 0, in n\ Uf =l B Lss! .( Zj ), 
we find that the constant bj t h, corresponding to Is + Rs(ou) satisfies 



As a result, 



M^llnc^W^ / 

i,h \j=l J Bl "iJ z ^ 



\\Qs(l s + Rs{co))\\ P 
<\\h + R8(oo)\\ p + cJ2\hh\ 

i,h 

<C\\l 5 \\ p + C\\Rs(u) 



dP s ,z,i 



dz 



i.h 



\k + R s (u)\ 



dz 3,h 



On the other hand, from Lemma [A. II and (|2.12p . we can deduce 



<5 \\p 



\P 

Zj,ag,i a S,j) + 



Cs s - 



3=1 

--0 



p-1 



p 



p-1 I 1 

\ln8\— + p 



For the estimate of H-R^o;)!^, we have 



m 

*«j { Ps > z ~ K i 



in 



2nq(x) 
I lnel 



m 
3=1 



2 r wq(x) 



3=1 



Kj Thi7 



<^TC||u; 

3=1 



\<>, Ps.Z ~ Kj - 



)y-i 






p 


27rg(x) N 


V 


| lne| / 





(3.16) 



Jln<J| 



P-3^ 



Thus, we obtain 
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G*(w)||oo<- { [\\k\\ P + 



\\RM\\ P ) 



.C\\d.5\—? 

()>' 

<r\ i ( 5 4 HI™ ^ (3 - 17) 



In 5 1 2 p \\nd\ 



Thus, G5 is a map from M to M. 

Step 2. G,5 is a contraction map. 

In fact, for any a;, G M, z = 1, 2, we have 

- G s u 2 = (QsLsy 1 Qs{Rs{^i) - Rgfa)). 

Noting that 

RsM = RsM = 0, in fi \ U J m =lJ B isa ,( Zj ), 
we can deduce as in Step 1 that 

(p-i) 2 
CI In 51 p 

||G?ia;i - GiWalU < 1 ll-Rtf(^i) - ^5(w 2 )|| P 

< c|k , r (JM^ + ^ |k _ W2|U 

<C5| hi ^^Hcdi - CJ 2 ||oo < - W 2 ||oo- 

Combining Step 1 and Step 2, we have proved that Gs is a contraction map from M to 
M. By the contraction mapping theorem, there is an unique u>s G M, such that us = Gsoj$. 
Moreover, it follows from (I3.17P that 



|w 4 ||qo < 5|ln<5| P 2 . 



4. Proof of The main results 



□ 



In this section, we will choose Z, such that Y^Jj=i Ps,zj+^s, where u$ is the map obtained 
in Proposition 13.31 is a solution of f ll.lOp . 
Define 

„ ■, P f,„ ,2 V- 1 [ ( 2m(x)Y +1 

m = - j a \m -Y,—J a x* (« - * - w ) + 

and 

K(Z) = I(P s>z + co s ). (4.1) 
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It is well known that if Z is a critical point of K(Z), then Y^jLi Ps,z,j + ^ is a solution of 

dnnj). 

In the following, we will prove that K(Z) has a critical point. 



Lemma 4.1. We have 



K{Z) = l[Y,Ps, 



5,Z,j 



o 



0=1 



(— 



Proof. Recall that 



We have 



K{Z) =l{P 5 ,z)+ [ 5 2 DP 5 ,zDlus + ^- [ \Du 5 \ 2 



m 

— Y 



Ps,Z + 0J5- Kj 



2nq(x) 



p+1 



2irq(x) 



p+i 



Using Proposition 13.31 and (I3.14p . we find 



P&.Z + — Kj 



2%q(x) 



In el 



p+i 



2nq(x) 



hid 



p+i 



-O 



ho 

'<5,i 



I lne|P 
On the other hand, 



P&.Z + — Kj 



O 



2nq(x] 
I In el 



p+i 



2nq(x) 
I In el 



p+i 



|lne|* 



<5 2 / DPs, z Dus = J2 / (W, 



/ (^^.ajj - a 5,i)+^5 

-3 

=0' 



\P 



|lne|* 
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Finally, we estimate 5 2 f n | -Da;^ | 2 . 
Note that 



+ U)$ — Kj — 



-5 2 Au s =J2xn 3 (Ps,z 
j=i ^ 

3=1 h=l V ° Z * h 



2irq(x 
'lne| 

's,z,j ' 



\ p m 

/ + j= i 



Hence, by ( 12.1 2p . we have 

5 2 [ \Du s \ 2 =jr [ (Ps, z + 
Jn . j Jo, [ V 



LUs — Kj 



27rq(x) 
I In el 



m 2 

+ EE 1 * 

m 



° z j,h 



3=1 



V|lne| 



m 2 



£ |foj,fel Halloo 

| lne|f 



--0 



So we can obtain that 



Lemma 4.2. VFe /iat>e 



dzL h dzi h 



i TPs. 



vi=i 



Proof. First, we have 



dK{Z) 
dz ijh 
d 



dPs,z <9w, 



dzi 



E 



27rg(z)V 



In el 



llnel^" 1 



(4 



ap 

<9z 



(5,Z 
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Since u$ G Eg )Z , we have 
Differentiating the above relation with respect to z^h, we can deduce 



r(p St z+u) S ).. 



du 5 



m 2 

EE'.* 



-5 2 A 



dP> 



S,Z,j 



dux 



dzjfi J dz i;h 

V>v^ , f fu r , P -i fdW 5<z a da s>j \ du s 



m 2 



|lne|P 



C— 



On the other hand, using (13.161) (for the definition of Rs(cu), see ( I3.12p ). we obtain 



E 

m 

E 



j=i J ^ 



27rg(a;) V / 27rg(ar^" 

n,z + ^ - — r-j — r - - /%• - 



|lne| 
27rg(a;] 

p-i 



dP> 



S,Z/, 



p 



- I -P<5,z _ S 



/ 2vrg(^)\ 



+ 

dP e ,z,< 

dz i>h — j Ul 



I In er| 

2vrg(x) x p 
I In el 



^5,Z — K j — 



dz 



i.h 



2ixq(x) 



p-1 



p-1 



(5,Z,-, 



'5J 



Inef 



Ps,z - K d - ) - \ W 5 , Zi ,a Sti - a 5A ) 



In el 



dP 



5.Z 



dzi 



(— 

\\kie\P 



O 



~~° ' \\ne\P~\ 
Thus, the estimate follows. 

Define 
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where rj > is a small constant and g > is a fixed small constant. Let 

D = {Z = Oi, ■■■ ,z m ) : Zi G fl e ,i = 1, • • • ,m, \zi - Zj\ > g L ,i^ j}, 

where fl g = [y : y G fl, d(y, dfl) > g\, and L > is a large constant. 
Denote K c = {Z : Z G D, K(Z) < c}. Consider 



dZ(t) 
dt 



-DK(Z(t)), t > 0, 



Z(0) G K c ^. 

Lemma 4.3. Z(t) does not leave D before it reaches K Cs > 1 . 
Proof. Note that 

, / \ 1 , 1 ft \ dh(x,z) 1 I x — z \ , t n . 

h(x,z) = —\n- q + ol, ^ ; =-— q (| -,n)+ol, 4.3 

27r F — z| an 27r|x — z| \ \x — z\ / 

if 2! is close to dfl, where n is the outward normal unit vector of d{x : x G fl, d(x, dfl) < 
<9fi)} and 2 is the reflection point of z with respect to dfl. 
Suppose that there is t > 0, such that Z(t ) =: (zi, ■ ■ ■ , z m ) 

(1) Suppose that there are i,j G {1, 2, ■ • • , m}, such that z 7^ j and — Zj| = g L '. 

Since cf(z, cftl) > £> and z ^ fl, using (14. 3p . we get < C'ln •= for any z, j, 



where C" > 0. Thus, we have 

G(zi,Zj) > In- — - — r -Cm- > Lin --Cm-. 

\Zi — Zj\ g g g 

Then, by Lemma 14.11 and Proposition \A.2\ we have 

C5 2 mC'5 2 In 1 L5 2 In 1 



if L > is large. 

(2) Suppose that there is i, such that Zi G <9fi e . 

Let n be the outward unit normal of dfl g at We have 

dG(zj,Zi) _ 1 / z^-^- \ dg(zj,Zi) 



dn \zj — Zi\ \ \zi — ZjY I dn 

where n is the outward normal unit vector of dfl g at Zi. 
On the other hand, if z,- G fi e , j 7^ z, satisfies 



0j Zj 



then, 
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/ *_SL n \ =o{\z i -z j \). 



So, we obtain 

/ 7.: — 7.z 

,n)> -C\zi - Zj\, V j ^i. 



Z{ Zj 



I Z i Z j I 



As a result, by Lemma [4.21 and Proposition IA.3[ we have 
OK A7r 2 5 2 Ki dqjzi) 2tt5 2 k 2 dg{zi, z { 

— i , ii, J? i I 



dn |lne||lnj| dn (hij) 2 dn 
^ 2n5 2 KiKj dg(zj, zj) C5 2 



(ln^) 2 dn (ln^) 2 ' 
On the other hand, we derive from ( 14. 3 p 

dg(zi,Zj) = 1 + o(l) 

and 

dg(zj,Zi) _ 1 + o(l) / Zj - Zj ^ 

where 2j is the reflection point of Z{ with respect to dfl. 

It is easy to check that if \zj — Zi\ < Mg, where M > is a fixed large constant, 
then 



I I 



dK > 2n6 2 ( K 2 + o(l) K^ + ojl) c >Q 



So 



dn ~ (In f ) 2 V 2g Mg 
Therefore, the flow does not leave D. □ 

Proof of Theorem \l.b\ We will prove that K(Z) has a critical point in K Cs - 2 \ K Cs < x . 

Suppose that K(Z) has no critical point in K Cs < 2 \ K 06 - 1 . Then from Lemma 14.31 that 
K 01 ' 1 is a deformation retract of K Cs < 2 . 

It is easy to see that K Cs < 2 = D and 

{Z : Z E D, \zi - Zj\ = g L , for some i ^ j] C if 04 - 1 . 
On the other hand, take i? large enough such that inf q > — V™ -, (g jjgi — V" 1 K »fffe> z j 
then K(Z) < implies that 
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(ln^) 2 ~ - |lne||ln^|' 



which implies that there are i ^ j, such that 



G(zj, Zi) > c In -, where d > is a constant. 



So, there is a a > 0, independent of 5, such that 

k - 2,- 1 < £ c 



Therefore 



: Z G D, \zi — Zj\ = g L , for some i ^ j} 
C K ** C {Z : Z E D,\ Zi - Zj \ < g a , for some i ^ j}. 

Filling the hole D* =: : Z E D , \zi — Zj\ = g L , for some i ^ j} in D, we obtain 



(4.4) 



{Z : Zj G \zi — zj\ < g L , for some i ^ j} 

C K^ 1 U D* C {Z : Zi G |2j - < for some i ^ j}. 

Since i*^ 1 is a deformation retract of K Cs < 2 , we find that il" c *4 U D* is a deformation 
retract of K Cs - 2 U D*. On the other hand, : z% G fi e , ^ = Zj, for some z 7^ j} is a 
deformation retract of : ^ G fi e , |zj — Zj| < g a , for some i 7^ j} if £ > is small. Using 
(14. 5p . we see that 

: Zi G = Zj, for some 2 7^ j} 



is a deformation retract of 



Q e x ••• x fi e = iT 5 > 2 UD*. 



This is impossible if f2 has nontrivial homology. 

i-p 

Thus we get a solution w$ for fll.lOp . Let u £ = ^^-ws, 5 = e yj^-j 2 , it is not difficult 

to check that u e has all the properties listed in Theorem II. 61 and thus the proof of Theorem 
11.61 is complete. 

□ 



Remark 4.4. In the proof of Theorem ll.6[ what we actually need is that the following 
function 



m m 



$(Z) = ^ 4n 2 Kiq(zi) + ^2-KK 2 i g(z i , z { ) -^irKjKjG^j, 

i=l i=l jjti 



REGULARIZATION OF POINT VORTICES 25 

as well as its small perturbation (in a suitable sense) has a critical point in D. Moreover, 
using the estimates as in Lemma 14.21 it is eas Y to see that if Y^j=iPs,z s ,j{ x ) + ^ is a 
solution of fll.lOp . and Z$ — > Zq as 5 — > 0, then Zq is a critical point of $>(Z). 



Proof of Theorem Note that the Kirchhoff-Routh function associated to the vortex 
dynamics is 



W(xi,--- ,X m ) = ^ KiKjG(Xi, Xj) + ^ ^ KjH(Xj, Xj) + y^Ktfl> (Xi). 

i^j i=l i=l 

Recall that h(zi, zj) = —H(zi, Zj), it is easy to check that 

rn 

$(Z) = -An 2 W(Z) + ™? ln R- 



i=l 



Hence, §(Z) and W(Z) possess the same critical points. 
By Lemma 14.11 14.21 and Proposition \A.2\ \A.3\ we have 



ln £ f-* 4(ln 

e i=l v e ' 

and 

dK(Z) 5 2 d$(Z) 

dzi h \lne\ 2 dzih 



+ 



d 2 


ln 


lne 




\ne 


3 



Thus, stable critical point of Kirchhoff-Routh function W(Z) implies that K(Z) has a 
critical point. So the result follows. □ 

Now we are in the position to prove Theorem 11.11 

Proof of Theorem \1.1[ By Theorem ll.4[ we obtain that u £ is a solution to (jl.9p . 
Set 

v e = (V(« 6 - g)) ± , u £ = V x v £ , 

then (v e , P £ ) forms a stationary solution for problem (11.11) . 
We now just need to verify 



» m 

/ uj £ — > 2_] Kj, as e — > 0. 
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By direct calculations, we find that 



Jci j=1 £ Jn 



Xn, [u e -q- 



Kj\ \ne\ \ p 
2ir 



= E 



\ne\ p f ( 2nq^ p 

W$ — Kj 



^(2tt)^7 y^ww "Mine,, , , 

JTt 111 

i=i S <S,3 
m 

— >■ ^ Kj, as e — > 0. 
i=i 

Therefore, the result follows. □ 

Remark 4.5. To regularize point vortices with equi-strength k, we do not need Xn P that 
is, we just need to consider the following problem 

'-e 2 Au = (u-q -£ln±£, x G Q, 
u — 0, a; G 90. 
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Appendix A. Energy expansion 

In this section we will give precise expansions of / (X]j=i P$,z,j^ an d af~-^ (Sj=i P&,z,^j > 
which have been used in section 4. 
We always assume that 

d(zj, d£l) > q > 0, \zi — Zj\ > g L , i ^ j, 
for some small g > and large L > 0. 

Lemma A.l. For x G % — 1, • • • , m, we have 

27iq(x) 

P8,Z(X) >Ki+ — , X G 5 a4i . ( i_ T84i .)(^), 

where T > is a large constant; while 
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27T(7 ( X ) 

Ps,z{x) < Ki + , x e Qi \ B SSi{l+sli) (zi), 

where a > is a small constant. 

Proof. Suppose that x G B Sgi (i-Ts Si )(zi)- It follows from (12.121) and <p[(s) < that 

2nq(x) ( ssi \ 

PsA x ) ~ K i~ .j. = W 5jZuas .{x) - a S]i + O I I 



|0'(l)||ln^fV s 5)i ; V|lne| 

if T > is large. On the other hand, if x G fij \ -B s ? (2A where a > o > is a fixed small 
constant, then 

_ , . 2nq(x) ^ . i? , . R 2'Kqix) 

P ^ (X) - Ki - "[hif = £ ^ ^ " " W + 

<Ca - ^ + o(l) < 0. 
Finally, if x G ^ s |.(^) \ £ S5>i (i+T S |.)(^) for some i, then 

27rg(x) / sJi 

- K i 1 ln I = W s>Zitasi {x) - a s>i + O I — -j- 



In 



= °*'* l n -Z- ~ as ' { + \ In JL 



ln(l + TsfJ / S f. 

if T > is large. Note that by the choice of a, B Sjj ^l+gj ^(^i) ^ -B Sil (i+Tsf f° r small 
5. We therefore derive our conclusion. □ 

Proposition A. 2. We have 

1 ^ S ' Z > ]„R + 2-^ Ah n R\2 + 2_^ I In c-l I In &\ fin ^2 



ln f tr 4 ( ln f) 2 tt iMiinfi tr ( ln f) 

nS 2 KiKjG(zj, Zi) f 5 2 ln I lne| 



where C is a positive constant. 
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Proof. Taking advantage of (12 . 3 p . we have 



Ja j=1 i=1 Jn 

m m „ / 

= J2J2 ( W s,w> - a sjY + W s , zuas . - rr^fo Zt 

j=l i=l J B aS]{Zj ) + \ In — 

First, we estimate 

/ (W StZi<a - a 5 ,i) P W s , Zi , as<i - 7^jr9{x, Zi 

{W s , Ziiasi - a 5y i) p+1 + a s> i / {W SlZi ,a S:i ~ as,i) P 

Bs Si {zi) JB as .(zi) 



(W 5tZit a Si - a Sti ) p g(x, 



r 2(p+l) „ r 2p /. 

(-) P l sl ^ +1 + a 5 J^-y- 1 slJ 

^Ss,i' ' JBi(0) ^S8,i' ^Bi(O) 



_ n(p+l) 5*aj ti grg^ 2*6^ ( £ 3 \ 
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By Lemma [O and f l27T2|) , 
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Thus, the result follows from Remark 12.21 



Proposition A. 3. We have 
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Proof. Direct computation yields that 
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Using (I2.12p , Lemma IA.1I and Remark 12. 2\ we find that 
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On the other hand, for Z 7^ i, from (12.121) . we have 



B Ls SJ (Zj) 



{Ws, Zj ,a s>j 



~ a S,jY + ~ ( P S,Z - Kj - 



2nq(x) 



8P; 



8,Z,l 



dz, h 



a S,j) 



p-1 Ss,j 



\\ns\ 



C 



In 



R 



--0 



|ln£|P+ 1 
Thus, the result follows. 
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